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I. INTRODUCTION 


Synthetic aperture radar (SAR) (i,2j is a side-looking sensor capable of 
very fine along-track resolution. This characteristic, in addition to its 
aii-weather monitoring capability, makes it an attractive instrument for 
generating space imagery. A SAR typically operates in a spectral region 
(1-10 GHz) that complements most optical scanners (LANDSAT). In this region, 
the subsurface sensing depth is much greater and the reflectance properties 
of imaged terrain objects are different [3]. A SAR will therefore collect a 
complementary set of "terrain signatures" that when registered with other 
spectral data can greatly facilitate the interpretation of remotely sensed 
imagery. But before multisensor data can be properly registered, the various 
types of distortion inherent in SAR imagery must be quantitatively analyzed. 
This requires an understanding of the properties of the sensor as well as the 
data processing system used in the image formation. 

The production of quality imagery from spaceborne SAR data requires 
extensive processing of the raw echo data. This procedure is more complex 
than previously used for aircraft SAR imagery due to the large increase in 
sensor altitude. The correlation procedure must compensate not only for 
undesirable spacec ift motion, but also for target motion resulting from 
Earth curvature and rotation [4,5]. To simplify this procedure, 
approximations are often made that can res ilt in both geometric and 
radiometric distortion in the image product. Furthermore, the viewing 
geometry of the SAR and the charac teristics of the data collection system 
will also introduce distortion into the imagery. Precise registration of SAK 
imagery with ether types of remotely sensed imagery requires the 
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identification of these distortions and the development of post--processing 
techniques to rectify them. Furthermore > considering the large quantity of 
data^ it is desirable that these po^<^^processing techniques be automated and 
designed to interface directly with the image processor to generate a 
geometrically and radi.ometrically correct product without supervision or 
operator interaction. 

This paper will summarize what is currently known about these 
distortions and describe the development to date of unsupervised 
post-processing rectification techniques* The geometric distortion can be 
divided into two categories. The first category consists of distortion 
derived from the radar viewing geometry. This includes such effects as 
ground range nonlinearities ^ radat foreshortening and radar layover. The 
second category consists of distortions introduced during the data 
processing* These distortions result from approximations made during the 
correlation such as in estimation of the target phase history, or 
compensation for the earth rotation* The processor induced distortions will 
obviously depend on the specific correlation algorithm used for image 
foraiStion* This paper generally addresses the effects on the image ^ roduct 
resulting from assumptions during the processing and it specifically 
considers distortions inherent in digital imagery produced by the digital 
image processor at JPL [6). 
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II. SAR REVIEW 


This section is intended as a brief tutorial on SAR system 
characteristics and processing techniques. Basic system parameters are 
introduced that are relevant to understanding the sensor and processor 

induced distortion to be discussed in subsequent sections. Interested 

# 

readers will find an excellent treatment of various characteristics of SAR in 
a recently published cwllection of papers 

2.1 SAR Sensor Characteristics 

The viewing geometry for a SAR in the imaging mode is illustrated in 
Figure 1. The sensor moves along a predetermined orbit track radiating 
energy in the form of pulses and receiving the backseat tered signal. This 
data is either recorded on film for later optical processing or is 
transmitted on a downiirk to a tracking station where it is digitized and 
recorded on high-density tape for digital processing. The antenna pointing 
direction is normal to the flight path and its orientation with respect to 
che spacecraft is norraally fixed. Two important system parameters are the 
antenna look angle 0 , defined as the angle the antenna beam with respect 
to nadir direction, and the slant range R, the distance from the sensor to the 
iraageo target area. The radar antenna range beamwidth (3^ and the 
spacecraft altitude H, determine the cross-track dimension of the radar 
footprint. 

The radar transmitter is typically designed to operate in a frequency 
region between L-band and X-band , thus permitting cloud and fog penetration. 
The pulse repetition frequency (FRF) of the radar and the spa .raft velocity 
determine the aziomth (along-t rack ) resolution. This is typically designed 
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to be several times higher than the range (cross-track) resolution so that 
the azimuth return can be divided into several looks for incoherent summation 
to reduce the speckle noise. The range resolution is determined by the 
spectral bandwidth c£ the transmitted pulse. Coded pulses ^uch as a chirp 
waveform are typically used to achieve the maximum resolution for a given 
level of transmitter power [8]. During the image processing these pulses 
must be compressed to a point response. Imrerfect compression can result in 
both reduced resolution and range sidelobes in the imagf . 

2.2 SAR Processor Characteristics 

Processing raw echo data into finished imagery essentially consists of 
the following four steps (16) • 

(1) Range pulse compression; 

(2) Doppler parameter estimation; 

(3) Azimuth correlation; and 

(4) Speckle reduction. 

2.2.1 Range Pulse Compression 

The purpose of this step is to compress the time dispersed radar 
transmitted pulse into an impulse. A chirp fanctioa is the radar pulse 
waveform typically used because of its focusing capability in an optical SAR 
processor. The key to attaining an adequate signal “to- e ratio is 
ti msmission of a wide“band pulse and use of the data processor to compress 
this pulse. The range resolution can be approximated by 

Pr ‘ 2W 
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where c is the speed of light and W is the pulse bandwidth* Spaceborne SAR 
requires unusually large pulse compression ratios compared with aircraft 
radar to achieve fine image quality* 

2*2*2 Doppler Parameter Estimation 

Formation of the synthetic aperture is critically dependent on accurate 
estimation of the phase delay history of the return echo signal. This phase 
delay can be approximated by (9) 


4)(t) 


^ 4 >( 0 ) + 




'd‘ * 2 


( 2 ) 


where 4>(0) is the phase at t=0, f^ is the Doppler frequency shift of the 

echv data* and f^ is the Doppler race, both evaluated at t=0* These two 

parameters, f . L.nd IV,, can be estimated based on the orbit ^.nd attitude 
d d 

data. If da accurate ephemeris is not available, these estimates car be used 
as initial predicts ana further refined by evaluation of certain 
characteristics in the resultant imagery. Details of these refinement 
techniques, termed autofocusing lor and clutterlock for can be 
found in the literature (9,10l. 

2.2.3 Azimuth Correlation 

One of the primary features of SAR is its very high solution, /z 
previously discussed, fine range resolution is achieved by transmission of a 
wide bandwidth pulse and compression of this pulse in the data processor. A 
high azimuth resolution is achieved by coherent processing of echo date from 
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successive r^dar pulses. The azimuth coire'.^uion operation essentially 
simulates a very narrow effective beamvidth in the azimuth direction by 
coherently adding the returns of several radar echoes. This operation is 
performed on range compressed data and inv^^.ves shifting the phase reference 
function given in (2) and summing th^ resultant detected echo data. 

The complexity of the azimuth correlation is increased by the range 
migration or th . target as it passes through the antenna beai-'. Tais 
migration consists of a range walk term which is directly proper r ror'il to th^t 
elapsed along-track time and a range curvature term, which is proportional to 
the square of the along-track time* To compensate for the that the 

target transverses may range resolution elements^ the inout data is resampled 
before * he phase reference function is appii» U. After phase detection, the 
return froji each pulse over the length of the synthetic aperture is summed to 
obtain a refined estimate of the target's brightness. If the returns are 
added over the fuxl aperture, a single-look full resolution image is 
produced. If the aperture is divided into sections, sevei^l single^looh 
reduced resolution images are produced. For example, twe single^look images 
can be produced by dividing the phase reference function in half and 
processing the tirst and second portions of the target response separately. 
The single-look images are then added incoherently as described in the next 
section to obtain a multiple-look xmage. 

2.2.4 Speckle Reduction 

Coherent processing of radar echo data is used to achieve very fine 
along-ti ck resolution, but this same proces.r»ing also makes th^- ’si^^ges 
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susceptible to speckling effects. Basically^ speckling results from 
scattering of notuuniformly distributed reflectors within a resolution cell 
during the period the cell is within the antenna beam. The speckles are 
signal dependent and therefore act like multiplicative noise. Techniques for 
speckle suppression essentially fall into two main categories: 1) averaging 

veral reduced resolution images produced from independent looks; and 
filtering a full-resolution single-look image to smooth the speckle. 

The multiple-look overlay technique is typically used for image 
production because it can achieve a satisfactory level of specicle reduction 
with a minimal amount of additional computation [11]. To produce multiple 
looks, the phase reference function is divided into segments during azimuth 
correlation. Each segment is applied to the appropriate portion of the 
return echo data to produce a reduced resolution single-look image of the 
target area. The single look imagery is then registered and averaged to 
produce a multi— look product with reduced speckle noise* 

The speckle can also be reduced by filtering a single-look high 
resolution image 112]. Smoothing algorithms have been devised based on the 
local statistics of the noise. In general, these filters can achieve a 
greater speckle reduction for a given resolution product then the 
multiple-look overlay technique, but because of the additional computational 
load, this approach is not feasible for large-scale image production. 
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111. GEOMETRIC DISTORTION 


SAR iinagery is subject to several types of geometric distortion that 
must be corrected before the imagery can be properly registered with other 
sensor data. These distortions can be classifed into two main categories: 

1) sensor derived distortion, and 2) processor derived distortion. 

3« i Sensor Derived Distortion 

The geometric distortions classified in this category are mainly derived 
from the viewing geometry of the radar. Design parameters such as the sensor 
altituue, look angle , and range beamwidth are critical in determining the 
degree of distortion in the imagery. In most cases, post-processing will be 
required to correct these distortion effects. 

3.1.1 Ground Range Nonlinearity 

Perhaps the most predominant distortion in a spaceborne SAR image is the 
ground range nonlinearity. This effect results from the fact that each pixel 
represents a constant slanc range distance rather than uniform ground 
spacing. Figure 2 illustrates the type of distortion inherent in a slant 
range projection. Features in the near range of a slant range image are 
compressed with respect to the far range. For example, a circle in the far 
range would have the appearance of an ellipsoid in the near range of the 
image. To project the linage in a ground range format, thtj ground distance 
represented by each pixel x.* the slant range projection must be determined 
Irom 


Ax 


A r 
sin4> 


(3) 
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where Ar is the slant range pixel spacing and ^ is the incidence angle of the 
radar beam at the target. The slant range spacing is constant and is given by 

Ar » ^ (4) 

s 

where c is the speed of light and is sampling frequency. The 
maximum frequency is limited by the range bandwidth. For SbASAT SAK, this 
frequency was 43.53 MHz producing a slant range spacing of 6.59 m. 

The incidence angle, 4> , can be related to the beam elevation angle, 9 , 
by trigonometric means as follows (Figure 3) 


4>* sin 


(H + R ) sin 0 
e 


R 

e 


(3) 


where is the radius of the earth, H is the spacecraft height and 6 , the 
look angle, is given by 


r 


9= cos 


-1 


+ (H + R ) R 
e e 

2R(H ^ R ) 
e 


( 6 ) 


The relationships for 4> and 0 in (5) and (6) have been derived assuming a 
smooth spherical surface with a radius equal to the radius of the earth at 
the target. Since the incidence angle is a function of the terrain, equation 
(3) is ony approximate for a rough surface. This could result in an error in 
the ground range projection. 
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3%1«2 Radar Foreshortening [13] 

Radar foreshortening is the variation in apparent size of identical 
features at different slopes with respect to the sensor. The effect for 
terrain sloping toward the radar is an increase in the effective incidence 
angle* Terrain sloping toward the radar will therefore appear elongated with 
respect to terrain sloping away from the sensor* This is illustrated in 
Figure 4. Since terrain features are recorded as a function of the slant 
range distance from the sensor, slope ab is mapped into a'b* in the image 
plane while be is mapped into b'c'* Although ab and be are identical in 
length, b*c* will be more than three times larger than a'b* in the image. 
Quantitatively stated the elongation (or shortening) factor is given by 


a*b* 

ab 


sin ( 6 • 3 ) 


17) 


b*c' 

be 


sin ( 0 ♦ i5 ) 


lb) 


where is the surface slope and 0 is the look angle* Thus for a look 
angle similar to the terrain slope, the foreshortening effect is most severe* 

The foreshortening can be corrected if the radar altitude and look angle 
are known and a terrain map of Che target area is available* One technique 
to correct the foreshortening is as follows (14): A digital topographica i 

map is illuminated from the same angle that the radar used to generate the 
original foreshortened image* This can be routinely done using a radar 
simulation program. The resultant movement of each pixel in the simulated 
image is recorded. This image is then registered to the actual radar image 
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and the inverse of this movement is applied* The result of this type of 
process is shown in Figure 5* Note that the presence of certain geological 
features becomes readily appare<^^ t* This technique, although capable of 
producing an irtage without foreshortening, is a tedious process which 
requires exact registration of the simulated image with the actual radar 
image* In addition, this technique requires a priori terrain information for 
the imaged area which may not be available* 

3.1*3 Radar Layover [13] 

Radar layover is a distortion inherent in all radar imaging of irregular 
terrain. It results from an image pixel being placed in its cross-track 
location based on its range from the sensor. If the target area contains 
features with steep slopes, the top of the feature can be at a closer range 
than the bottom* This effect is shown in Figure 6* It is especially severe 
when features of appreciable level relief appear in the near range* 

3.1.4 Earth Rotation 

The earth rotation effect is common to satellite scanning systems where 
the scan lines are gathered serially in time while the earth's SUPfdCe IS 
rotating during the imaging period* This results in the near edge of the 
image being skewed with respect to the nadir track of the satellite as shown 
in Figure 7a. This effect, known as range walk or range migration can be 
compensated by resampling the echo data during the image processing. 

Figure 7b shows the outline of an image frame after data resampling to 
compensate for the range walk* The stairstep edge is an artifact of the 
muitiple-look overlay process and does not result in any discontinuities 
within the frame* 
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3.2 Processor Derived Distortion 


The distortions generated during image formation are obviously dependent 
on the specific processing algorithm* Generally » approximations made in 
determining the phase reference function will result in a skewed image* This 
section discusses some typical approximations and the resultant distortions* 


3*2*1 Earth Curvature and Rotation 

In the previous section, the range walk effect that results from the 

Earth rotation was described* This rotation in conjunction with the £irth 

curvature can also result in an image skew in the azimuth direction* This 

skew occurs during the image formation process* As previously discussed, to 

process raw echo data into imagery the Doppler frequency and frequency rate 

must be determined. The rotation and curvature of the earth, however, result 

in a complex Doppler response that is a function not only of target latitude 

but also depends on target position within the swath* A plot of Doppler 

frequency as a function of swath position and latitude is s^ wn in Figure 8. 

It is extremely difficult to design a processor that adapts to the varying 

frequency without introducing reduced image resolution resulting from 

misregistration during the multiple^look overlay* As an alternative to 

precise tracking of the Doppler frequency during the processing, an optimum 

f is selected for the entire image frame and ^he resultant skew is 
d 

corrected in the post-processing* This skew resulting from the Doppler 
mismatch is essentially in the azimuth direction as shown in Figure 9. Th^^ 
amount of skew is dependent on the difference between the actual Doppler 
response of each target and the f^ used in the processing and is given by 


AN 


at 



C9) 
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where PRF is the pulse repetition frequency, L is the number of looks, Af^ 

is the Doppler frequency mismatch, f . is the Doppler rate and AN is the 

azimuth displacement in pixels* Note that since the iso-^Doppler lines arf^ 

nearly linear» a linear correction can be used to deskew the image* For a 

lOC^ km swath width as in SEASAT SAR, a typical skew is AN 150 pixels or 

az 

2*5 km* which corresponds to a skew angle of 1*5 degrees* 

3*2*2 Sensor Parameter Shift 

The eccentricity of the spacecraft orbit requires adjustment of the 
sensor parameters at certain intervals during the orbit* This can result in 
either a discontinuity} or variation of the pixel resolution within an image 
frame. A sharp discontinuity in the image results from a shift in the delay 
of the pulse sampling window* This parameter shift is made in the data 
collection system to maintain the optimum signal* ter noise ratio in the 
sampled echo data* If the raise sampling window were not adjusted as the 
spacecraft altitude changed, this window would not center on the portion of 
the signal with the strongest return and the resulting imagery would be 
degraded* The effect on an image produced from data collected during a 
sampling window shift is to displace one part of the image in range with 
respect to the other, resulting in a line of discontinuity across the frame* 

Another parameter requiring adjustment based on variation in the 
spacecraft orbit is the pulse repetition frequency (PRF). The PRF is 
adjusted to prevent pulse mixing or range ambiguities* An upper limit on the 
PRF is set to avoid interference between successive pulses* This limit is 
dependent on the range swath width which is in turn a function of spacecraft 
altitude* A lower limit on the PRF is determined by the resolution and 
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azimuth ambiguity requirements. This bound is also dependent on the 
spacecraft orbit. To keep the PRF within these bounds^ it must be adjusted 
according to variation in the spacecraft height. Certain frames therefore 
may contain pixels of mixed resolution if the PRF was adjusted during that 
frame's data collection period. 

3.3 Radiometric Characteristics 

Registration of a SAR image to an optically sensed image or another SAR 
image can present a difficult problem since the apparent brightness of an 
image pixel is critically dependent on the relative position of the sensor. 

A small change in the aspect angle can significantly alter the appearance of 
a scene as shown in Figure 10. Furthermore, SAR imagery is corrupted by both 
speckle and thermal noise, pulse compression sidelobes> ambiguity responses 
and weak signal suppression effects [6]. With proper radar and processor 
design, the effects of most of these distortions on the final image product 
can be minimized. The image degradation resulting from speckle noise is 
greatest in coherent imagery of terrain with a surface roughness comparable 
to the wavelength of the illumination [13]. The speckle can be reduced by 
summation of looks or filtering at the cost of pixel resolution. 

When attempting to register SAR imagery from two adjacent passes for the 

purpose of constructing a mosaic, the distortion resulting from speckle can 

cause misregistration. This results from the fact that for a small change in 

the sensor ^^osition (4 km for SEASAT) the speckle noise in the resultant 

imagery is totally independent. Therefore, from one pass to the next a 

feature may be distorted differently resulting in a poor correlation between 

the two frames. An additional factor that can cause image misregistration 

between two adjacent passes is a phenomenon called specular point 
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migration illustrated in Figure ll« This effect occurs wh^jn a feature with a 
varying slope is imaged from two different sensor positions » The apparent 
position of a feature changes relative to other features within the frame 
because the point of maximum reflectivity of that feature occurs at a slope 
perpendicular to the radar beam. This effect can also result in apparent 
misregistration of features in SAR imagery. 
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IV* CURRENT ?KOGhESS IN REGISTRATION/RECTIFICATION OF SAR FREQUENCY 


Many of the geometric and radiometric distortions present in digitally 
processed SAR imagery do not occur in optically sensed imagery* The effect 
of these distortions on image resolution, rectification and ultimately 
mosaicking have not been fully investigated* To date very few mosaics have 
been constructed with digital SAR imagery and at present no comprehensive 
studies have been conducted to investigate the feasibility of generating 
unsupervised mosaics* 

The SAR processing research group at Jet Propulsion Laboratory has laid 
the ground work for development of an unsupervised mosaicking procedure and 
plans to implement this procedure within the next year* Recently completed 
is an algorithm capable of determining the absolute location of an image 
pixel without the aid of ground reference points [6]* This algorithm 
utilizes information provided by the spacecraft ephemeris and characteristics 
of the sensor data collection system to predict the latitude and longitude of 
an arbitrary image pixel* Tests have chovn this technique has an accuracy of 
better than 200 m for SEASAT SAR data. The target location uncertainty is 
primarily dependent on the accuracy of the spacecraft position and velocity 
provided by the ephemeris and the validity of the assumed geoid in the target 
area* 


The capability for unsupervised pixel location is a prerequisite for 
automated rectification and mosaicking of SAR imagery* With informtion on 
both Che sensor and target position, the parameters for slant range to ground- 
range conversion can be determined as well as the Doppler shift in the echo 
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data. Using this information and the image processing parameters, it is a 
simple procedure to conduct range and azimuth interpolation to remove the 
ground range nonlinearity and azimuth skew distortions* An algorithm has 
been developed at JPL to correct these distortions for digitally cor-- ‘ d 
SEASAT imagery. A preliminary version of this algorithm has been piemen’^ 
utilizing the hardware architecture of the digital processor. This hardware 
configuration, shown in Figure 12, was designed both for efficiency and 
economy [17]. The host computer is a Systems Engineering Laboratory 32/77 
which features very rapid I/O data transfers. Attached to the host are three 
Floating Point Systems AP120b array processors and four 300 Mbyte discs. The 
image data is transferred from disc to the array processors through the 
host. Each AP operates simultaneously on a separate area of the image and 
returns the rectified product to the host to be merged and transferred to 
disc. It is anticipated that with proper design this approach can rectify a 
6000 X 6000 pixel frame which covers a 100 x 100 km ground area in under 
30 minutes. This approach is compatible with the image processor and could 
be easily incorporated into the image production procedure. 

A prelir^inary version of this algorithm is curr* vtly operating at JPL. 

It utilizes a single AP and makes some simplifying assumptxons in generating 
the rectification parameters. The current operating time for a full frame 
is 80 minutes. The output has been resampled to 12.3 m spacing in both the 
range and azimuth dimensions. As a result of approximation, the rectified 
product still contains some residual skew. Tests show approximately a 200 ra 
registration error across a lOO'-km frame. The expected registration error 
using a more exact version of the algorithm is less than 50 m or about 
2 pixels. 
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Following completion and test of the rectification algorithm, studies 
are planned to investigate the feasibility of generating unsupervised image 
mosaics. Imagery from parallel passes as well as ascending and descending 
passes will be examined. Statistics will be compiled on the correlation of 
common features as well as the misregistration between geometrically 
corrected images from two adjacent passes. Automated feature detection 
techniques are also being considered to locate suitable tie points. 

Obviously much work remains to be done with regard to rectification and 
registration of SAR imagery. It is clear that digitally correlated imagery 
offers distinct advantages over optical imagery such as: reproducible 

processing, large inherent dynamic range and highly accurate pixel location. 
The refinement of these techniques utilizing SEASAT SAR data will provide 
insight into the design of future sensor systems to optimize the registration 
and rectification accuracy. 
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Figure 1. Viewing geometry for a typical side-looking synthetic, aperture, radar 
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Figure 2. Slant range to ground range nonlinearity inherent 
in all side looking radat systems. Near range 
features are compressed relative to far range 
features because of change in incidence angle, 4i. 


187 



SENSOR 


OK 

OF 


POOR QUALtTY 



Fik;ure 3. Rolationshlp between look an^le, . and incidence 
angle c tor smooth spherical planet surtacc. 
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Figure 5. Kesulls ot ^ (>F'ipop.s;it i on i'^M* r.ular t'oroshor te^^.ing using 
digital terrain mode 1 . :\) Simulated radar image from 

t - ' p o g r a p h i e. a I ma p , b) C v ^ r r e e. L ed. S I' A S A i ima ge , c ) Origi- 
nal SKASAT image. NvUe the erdnatu' emeu t ol taull lines 
not evident in original image toviiaghi^ 1^81). 
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Figure 6, Radar layover is a function uf the terrain 

slope and the slant range to the target. Layover 
occurs when the surface feature slope exceeds the 
look angle. Signal scattered from the top of 
features 1 and 2 vill be placed at a nearer cross- 
track position in the image than signal scattered 
from the bottom. 
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a 
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Figure 7. Range migration of the ground track resulting from 
earth rotation during the imaging period is shown 
in a. The range skew distance is given by d^ = ^^c***^ 
W(K-‘ro is the velocity of the earth and At is the 
imaging time for the frame. The cross track skew 
shown in b is applied during the image formation. 



DOPPLER FREQUENCY, Hz 



SLANT RANGE, km 


Figure 8. Plot of Doppler frequency as a function of 

swath position and latitude • Nominal values 
for SEASAT SAR were assumed for spacecraft 
in descending mode (Wu, 1981) • 
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Figure 10. Effect of radar aspect angle on the apparent 
brightness of target. Figure ICa was imaged 
during a descending mode and Figure 10b during 
an ascending mode 
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Figure 11* Specular point migration results when sensor images 
the same feature from two different positions. 

Point of maximum return is dependent on target slope 
and may vary relative to other features within the 
fra a. 
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Figure 12. Block diagram of digital processing system at 
Jet Propulsion Laboratory. Data is recorded 
on high density digital tape (HDDT) at the 
tracking station. It is received at the 
processing facility via an optical link and 
recorded on disk for processing (Wu, 1981) . 
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